Improved detection of artifactual viral minority variants in high-throughput sequencing data

High-throughput sequencing (HTS) of viral samples provides important information on the presence of viral minority variants. However, detection and accurate quantification is limited by the capacity to distinguish biological from artificial variation. In this study, errors related to the Illumina Hiseq2000 library generation and HTS process were investigated by determining minority variant frequencies in an influenza A/WSN/1933(H1N1) virus reversegenetics plasmid pool. Errors related to amplification and sequencing were determined using the same plasmid pool, by generation of infectious virus using reverse genetics followed by in duplo reverse-transcriptase PCR (RT-PCR) amplification and HTS in the same sequence run. Results showed that after 'best practice' quality control (QC), within the plasmid pool, 1 minority variant with a frequency >0.5% was identified, while 84 and 139 were identified in the RT-PCR amplified samples, indicating RT-PCR amplification artificially increased variation. Detailed analysis showed that artifactual minority variants could be identified by two major technical characteristics: their predominant presence in a single read orientation and uneven distribution of mismatches over the length of the reads. We demonstrate that by addi

Improved detection of artifactual viral minority variants in high-throughput sequencing data

High-throughput sequencing (HTS) of viral samples provides important information on the presence of viral minority variants. However, detection and accurate quantification is limited by the capacity to distinguish biological from artificial variation. In this study, errors related to the Illumina HiSeq2000 library generation and HTS process were investigated by determining minority variant frequencies in an influenza A/WSN/1933(H1N1) virus reverse-genetics plasmid pool. Errors related to amplification and sequencing were determined using the same plasmid pool, by generation of infectious virus using reverse genetics followed by in duplo reverse-transcriptase PCR (RT-PCR) amplification and HTS in the same sequence run. Results showed that after "best practice" quality control (QC), within the plasmid pool, one minority variant with a frequency >0.5% was identified, while 84 and 139 were identified in the RT-PCR amplified samples, indicating RT-PCR amplification artificially increased variation. Detailed analysis showed that artifactual minority variants could be identified by two major technical characteristics: their predominant presence in a single read orientation and uneven distribution of mismatches over the length of the reads. We demonstrate that by addition of two QC steps 95% of the artifactual minority variants could be identified. When our analysis approach was applied to three clinical samples 68% of the initially identified minority variants were identified as artifacts. Our study clearly demonstrated that, without additional QC steps, overestimation of viral minority variants is very likely to occur, m

Companion Animals as a Source of Viruses for Human Beings and Food Production Animals

Companion animals comprise a wide variety of species, including dogs, cats, horses, ferrets, guinea pigs, reptiles, birds and ornamental fish, as well as food production animal species, such as domestic pigs, kept as companion animals. Despite their prominent place in human society, little is known about the role of companion animals as sources of viruses for people and food production animals. Therefore, we reviewed the literature for accounts of infections of companion animals by zoonotic viruses and viruses of food production animals, and prioritized these viruses in terms of human health and economic importance. In total, 138 virus species reportedly capable of infecting companion animals were of concern for human and food production animal health: 59 of these viruses were infectious for human beings, 135 were infectious for food production mammals and birds, and 22 were infectious for food production fishes. Viruses of highest concern for human health included hantaviruses, Tahyna virus, rabies virus, West Nile virus, tick-borne encephalitis virus, Crimean-Congo haemorrhagic fever virus, Aichi virus, European bat lyssavirus, hepatitis E virus, cowpox virus, G5 rotavirus, influenza A virus and lymphocytic choriomeningitis virus. Viruses of highest concern for food production mammals and birds included bluetongue virus, African swine fever virus, foot-and-mouth disease virus, lumpy skin disease virus, Rift Valley fever virus, porcine circovirus, classical swine fever virus, equine herpesvirus 9, peste des petits ruminants virus and equine infectious anaemia virus. Viruses of highest concern for food production fishes included cyprinid herpesvirus 3 (koi herpesvirus), viral haemorrhagic septicaemia virus and infectious pancreatic necrosis virus. Of particular concern as sources of zoonotic or food production animal viruses were domestic carnivores, rodents and food production animals kept as companion animals. The current list of viruses provides an objective basis for more in-depth analysis of the risk of companion animals as sources of viruses for human and food production animal health

The Carriage of Multiresistant Bacteria after Travel (COMBAT) prospective cohort study: Methodology and design

Background: Antimicrobial resistance (AMR) is one of the major threats to public health around the world. Besides the intense use and misuse of antimicrobial agents as the major force behind the increase in antimicrobial resistance, the e

Genetic and antigenic characterization of influenza A/H5N1 viruses isolated from patients in Indonesia, 2008–2015

Since the initial detection in 2003, Indonesia has reported 200 human cases of highly pathogenic avian influenza H5N1 (HPAI H5N1), associated with an exceptionally high case fatality rate (84%) compared to other geographical regions affected by other genetic clades of the virus. However, there is limited information on the genetic diversity of HPAI H5N1 viruses, especially those isolated from humans in Indonesia. In this study, the genetic and antigenic characteristics of 35 HPAI H5N1 viruses isolated from humans were analyzed. Full genome sequences were analyzed for the presence of substitutions in the receptor binding site, and p

Prevalence and risk factors for carriage of ESBL-producing Enterobacteriaceae in a population of Dutch travellers: A cross-sectional study

Background: We investigated prevalence and predictive factors for ESBL-E carriage in a population of mostly travellers prior to their travel (n = 2216). In addition, we examined ESBL genotype before travel and compared these to returning travellers. Method: A questionnaire and faecal sample were collected before travel, and a second faecal sample was collected immediately after travel. Faecal samples were analysed for ESBL-E, with genotypic characterization by PCR and sequencing. Risk factors for ESBL-E carriage prior to travel were identified by logistic regression analyses. Results: Before travel, 136 participants (6.1%) were colonized with ESBL-E. Antibiotic use in the past three months (ORadjusted 2.57; 95% CI 1.59–4.16) and travel outside of Europe in the past year (1.92, 1.28–2.87) were risk factors for ESBL-E colonisation prior to travel. Travel outside of Europe carried the largest attributable risk (39.8%). Prior to travel 31.3% (40/128) of participants carried blaCTX-M 15 and 21.9% (28/128) blaCTX-M 14/18. In returning travellers 633 acquired ESBL-E of who 53.4% (338/633) acquired blaCTX-M 15 and 17.7% (112/633) blaCTX-M 14/18. Conclusion: In our population of Dutch travellers we found a pre-travel ESBL-E prevalence of 6.1%. Prior to travel, previous antibiotic use and travel outside of Europe were the strongest independent predictors

Carriage of Blastocystis spp. in travellers - A prospective longitudinal study

Introduction: A lack of prospective and longitudinal data on pre- and post-travel carriage of Blastocystis spp. complicates interpretation of a positive test post-travel. Therefore we studied dynamics of Blastocystis carriage in a cohort of Dutch travellers. Methods: From the prospective, multicentre COMBAT study among 2001 Dutch travellers, a subset of 491 travellers was selected based on travel destination to 7 subregions (70 or 71 travellers each). Faecal samples taken directly before and after travel were screened for Blastocystis with qPCR, followed, when positive, by sequence analysis to determine subtypes. Results: After exclusion of 12 samples with missing samples or inhibited qPCR-reactions, stool samples of 479 travellers were analysed. Before travel, 174 of them (36.3%) carried Blastocystis and in most of these, the same subtype was persistently carried. However, in 48/174 of those travellers (27.6%; CI95 20.8–36.6%) no Blastocystis or a different subtype was detected in the post-travel sample, indicating loss of Blastocystis during travel. Only 26 (5.4%; CI95 3.7%–8.0%) of all travellers acquired Blastocystis, including two individuals that were already positive for Blastocystis before travel but acquired a different subtype during travel. Discussion: This study shows that Blastocystis carriage in travellers is highly dynamic. The observed acquisition and loss of Blastocystis could either be travel-related or reflect the natural course of Blastocystis carriage. We demonstrate that the majority of Blastocystis detected in post-travel samples were already carried before travel

Global phylogenetic analysis of Escherichia coli and plasmids carrying the mcr-1 gene indicates bacterial diversity but plasmid restriction

To understand the dynamics behind the worldwide spread of the mcr-1 gene, we determined the population structure of Escherichia coli and of mobile genetic elements (MGEs) carrying the mcr-1 gene. After a systematic review of the literature we included 65 E. coli whole genome sequences (WGS), adding 6 recently sequenced travel related isolates, and 312 MLST profiles. We included 219 MGEs described in 7 Enterobacteriaceae species isolated from human, animal and environmental samples. Despite a high overall diversity, 2 lineages were observed in the E. coli population that may function as reservoirs of the mcr-1 gene, the largest of which was linked to ST10, a sequence type known for its ubiquity in human faecal samples and in food samples. No genotypic clustering by geographical origin or isolation source was observed. Amongst a total of 13 plasmid incompatibility types, the IncI2, IncX4 and IncHI2 plasmids accounted for more than 90% of MGEs carrying the mcr-1 gene. We observed significant geographical clustering with regional spread of IncHI2 plasmids in Europe and IncI2 in Asia. These findings point towards promiscuous spread of the mcr-1 gene by efficient horizontal gene transfer dominated by a limited number of plasmid incompatibility types

Informing epidemic (research) responses in a timely fashion by knowledge management - a Zika virus use case

The response of pathophysiological research to emerging epidemics often occurs after the epidemic and, as a consequence, has little to no impact on improving patient outcomes or on developing high-quality evidence to inform clinical management strategies during the epidemic. Rapid and informed guidance of epidemic (research) responses to severe infectious disease outbreaks requires quick compilation and integration of existing pathophysiological knowledge. As a case study we chose the Zika virus (ZIKV) outbreak that started in 2015 to develop a proof-of-concept knowledge repository. To extract data from available sources and build a computationally tractable and comprehensive molecular interaction map we applied generic knowledge management software for literature mining, expert knowledge curation, data integration, reporting and visualization. A multi-disciplinary team of experts, including clinicians, virologists, bioinformaticians and knowledge management specialists, followed a pre-defined workflow for rapid integration and evaluation of available evidence. While conventional approaches usually require months to comb through the existing literature, the initial ZIKV KnowledgeBase (ZIKA KB) was completed within a few weeks. Recently we updated the ZIKA KB with additional curated data from the large amount of literature published since 2016 and made it publicly available through a web interface together with a step-by-step guide to ensure reproducibility of the described use case. In addition, a detailed online user manual is provided to enable the ZIKV research community to generate hypotheses, share knowledge, identify knowledge gaps, and interactively explore and interpret data. A workflow for rapid response during outbreaks was generated, validated and refined and is also made available. The process described here can be used for timely structuring of pathophysiological knowledge for future threats. The resulting structured biological knowledge is a helpful tool for computational data analysis and generation of predictive models and opens new avenues for infectious disease research. ZIKV Knowledgebase is available at www.zikaknowledgebase.eu